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Using the AdS/CFT correspondence for strongly coupled gauge theories, we calculate thermaliza- 
tion of mesons caused by a time-dependent change of a baryon number chemical potential. On the 
gravity side, the thermalization corresponds to a horizon formation on the probe flavor brane in the 
AdS throat. Since heavy ion collisions are locally approximated by a sudden change of the baryon 
number chemical potential, we discuss implication of our results to RHIC and LHC experiments, 
to find a rough estimate of rather rapid thermalization time-scale tth < 1 [fm/c]. We also discuss 
universality of our analysis against varying gauge theories. 



Introduction 

The AdS/CFT correspondence [IH3], or more broadly, 
the gauge/gravity duality, is an extremely useful tool to 
study strongly coupled field theories. Recently, this cor- 
respondence has been applied to various field theory set- 
tings, and these applications open up many new corre- 
spondences between gravity to other branches of physics. 
Perhaps one of the most surprising things in the success 
of using gravity to study strongly coupled gauge theo- 
ries is that it seems to work even for an explanation 
of heavy-ion collision experiment data at Brookhaven's 
Relativistic Heavy Ion Collider (RHIC). In RHIC exper- 
iments |H [5] , one of the big surprises was that a quark- 
gluon plasma (QGP) forms at a very early stage [5] just 
after the heavy ion collision, i.e. a rapid thermalization 
occurs. This obviously requires a theoretical explana- 
tion, but remains as a challenge, because this requires a 
calculation of the strongly coupled field theory in non- 
equilibrium process. In this letter, we study the ther- 
malization in strongly coupled field theories by using the 
gauge/gravity duality. 

The key idea is to approximate the heavy ion col- 
lision by a sudden change of a baryon-number chemi- 
cal potential locally at the collision point. Using the 
AdS/CFT correspondence, we obtain strongly coupled 
gauge theory calculations for the thermalization, where 
a time-dependent confinement /deconfinement transition 
occurs due to a sudden change of the baryon-number 
chemical potential, with dynamical degrees of freedom 
changing from mesons to quark/gluon thermal plasma. 
We calculate a time-scale for that. Our strategy can 
be summarized briefly as follows; On the gravity side 
of the AdS/CFT correspondence, the change in the 
baryon chemical potential is encoded in how we throw in 
the baryonically-charged fundamental strings (F-strings) 
from the boundary to the bulk. Since the F-string end- 
point is a source term for the gauge fields on the flavor 
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brane in the AdS bulk, this provides a time-dependent 
gauge field configuration. This induces a time-dependent 
effective metric for the degrees of freedom on the flavor 
brane, which are mesons. As a result, this yields the 
emergence of an apparent horizon on the flavor brane, 
which signals, in the dual strongly coupled field theory, 
the "thermalization of mesons" , which we mean that the 
meson degrees of freedom change into quark and gluon 
degrees of freedom with thermal equilibrium. 

Any computation of thermalization of mesons due to 
the injection of the baryon charge in strongly coupled 
gauge theories has never been proposed. We provide a 
generic framework for it in this paper. We also present 
computations at different gauge theories, and argue how 
universal the thermalization time-scale is. We also dis- 
cuss its implication in strongly coupled gauge theories, 
which hopefully offers a path to realistic QCD. The obser- 
vation of the flavor thermalization due to changes of ex- 
ternal parameters (i.e., quantum quench) different from 
the baryon charge, was studied in [J]. Previous stud- 
ies on holographic thermalization, for example in [8Hl2j. 
discussed glueball sectors, while ours observes the meson 
sector thermalization. What we see is the thermaliza- 
tion on the probe flavor brane. Since the back-reaction 
is not taken in our setting, this thermalization is not at 
all related to the one of the glueball sectors. 

In our framework, the only input is the function which 
represents how we throw in the baryonically charged F- 
strings. Therefore, we have small number of parameters, 
which includes a typical maximum value of the baryon 
density and the time-scale for changing the chemical po- 
tential. With collision parameters at RHIC, we obtain 
the thermalization time-scale as t t h < 1 [fm/c]. Actually 
this time-scale can be well compared with the known hy- 
drodynamic simulation requirement i t h < 2 [fm/c] dis- 
cussed for example in (T3h17| . We also "predict" that 
heavy-ion collisions at CERN's Large Hadron Collider 
(LHC) exhibit slightly smaller order of the time-scale for 
thermalization as i t h ^ O (0.1) [fm/c]. 

Let us make a few more comments about comparing 
our analysis with data. As we mention previously, we are 
discussing the time-scale of "thermalization of mesons" , 
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which is the horizon formation on the probe brane only. 
In the real world experiments such as RHIC or LHC, 
thermalization should involve not only mesons but also 
glueballs, which corresponds to the black hole horizon 
formation in the bulk geometry, not only on the probe 
brane. Since we are not treating the glueball sectors, the 
"prediction" above in order to compare our analysis with 
the data is unfortunately not as accurate as it should be. 
Another worry for comparison with data is that we are 
discussing the large TV limit of gauge theories. Therefore, 
the reader should regard our analysis as just an indication 
of rapid thermalization of some sectors of the large N 
gauge theories. 

In the following, after solving the equations in the 
gravity side with a generic time-dependent baryon chem- 
ical potential, we compute the apparent horizon and the 
time-scale for the thermalization. The simplest example 
offered is Af = 4 super Yang-Mills with Af — 2 hypermul- 
tiplets as "quarks". We conclude with the statement of 
the universality, by showing some variations of the setup, 
including quark masses and confining scales. 

D3-D7 system with quark injection 

The simplest set-up in AdS/CFT with quarks is the 
Af = 2 supersymmetric massless QCD constructed by a 
D3-D7 system [TS], where we consider the gravity back- 
reaction of only D3-branes and regard D7-brane as a 
prove flavor brane. We are interested in the dynamics of 
mesons and the deconfinement of quarks, which is totally 
encoded in the probe flavor D7-branc in the AdS$ x S 5 
geometry, 
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FIG. 1: How the baryon charge density changes in time, as we 
throw-in fundamental strings from the boundary. Top figure: 
there is no baryon charge. Bottom figure: the static baryon 
(quark) charge is provided by the fundamental strings (dashed 
lines) connecting the probe D7-brane and the AdS horizon 
(black blob). The end points of the fundamental strings are 
electric charges (flavor charges) on the D7-brane. An electric 
flux (solid lines with arrows) emanates from those charges 
on the D7-brane. Middle figure: the flavor electric charges 
are thrown-in from the AdS boundary into the bulk. This 
specifies the time-dependent chemical potential. The electric 
flux is time-dependent. The big arrow denotes the motion of 
the electric charges on the D7-brane. 
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is the AdS radius defined by R 4 = Airg s N c a' 2 . The string 
coupling is related to the QCD coupling as 2itg s = <?q CD • 
The dynamics of the flavor D7-brane is determined by the 
D7-brane action 



S = -/i 7 J d 8 ^-dct(G ab + 2na'F ab 



(2) 



where the D7-brane is at w§ — and are extended on the 
gauge theory directions x^ (// = 0, 1, 2, 3), p, and fi 3 . The 
fluctuations of gauge fields A a (or, scalar field -q = wq) on 
the D7-brane corresponds to vector (or, scalar) mesons. 
For a concise review of the D3-D7 system and meson 
dynamics, see [19 . G a b is the induced metric on the 
D7-brane. The asymptotic (p — > 00) value of r\ corre- 
sponds to the quark mass m q — r\ / '2ti a' . For simplicity 
we first put it zero. This corresponding to a "marginal 
confinement" for the mesons on the flavor brane since it 
has only zero-sized horizon [27) . The D7-brane tension is 
p 7 ee 1/(2tt)V*' 4 - 

In AdS/CFT, the response to the change in the baryon 
chemical potential is totally encoded in this D7-brane 



action. We will solve this gauge field for arbitrary time- 
dependent chemical potential. 

As the chemical potential of our interest is homoge- 
neous, we may turn on only the F tr component. With a 
redefinition of the AdS radial coordinate z = R 2 /p, the 
D7-brane action ^ is equivalent to a 1+1-dimensional 
Born-Infcld system in a curved background, 



S = -/J 7 ^Vol(S* 3 ) Jdtdz 



1-— (2™')^ 



2 

tz ' 



(3) 



where V3 is the volume of x 1 ' 2 ' 3 space. 

In the AdS / CFT dictionary, the static baryon chemical 
potential corresponds to A t (r — 00) — A t (r = 0). This 
basically counts the number of electric charges located 
at the origin r = 0. In order to change this number in 
a time-dependent manner, we need to consider an addi- 
tional source term 



SS = p 7 V 3 Vo\(S 3 ) J dtdz (A t f + A, 



J" 



(4) 



which describes end-points of a fundamental string (elec- 
tric charges) thrown in from the outside of the system, 
i.e., from the boundary into the bulk. See Fig. [I] Since 
the geodesic of the fundamental string end-points is the 
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light geodesies determined by the induced metric G a b, 
and in the present case Gu — —G zz , it is just along a null 
vector (v t ,v z ) = (1, —1). Therefore the source current is 
an arbitrary function of the variable t — z. With a cur- 
rent conservation relation, we obtain that j = j z , and we 
take the arbitrary source function as j* = j z = g'(t — z). 
Given this source current j, the gauge field strength F tz 
is readily solved from the equations of motion on the D7- 
brane and we obtain 



(2ira r )F tz = 



R 2 z g(t-z) 



^/(27ra') 2 i? 12 + z e (g(t-z)Y 



(5) 



This is the gauge field solution which encodes the infor- 
mation of the time-dependent chemical potential given 
the source current j* = j z = g'(t — z). 

The relation between g and the chemical potential is 
linear. For the static case g(t—z) = constant, the solution 
([5| is nothing but a conventional Born-Infeld solution on 
a D-brane. So we can compare the Born-lnfcld charge 
g with the fundamental string charge (which equals the 
quark number), following the techniques found in [20], to 
obtain 



g(t) = (2/n)(2Tra r ) 4 \n B (t). 



(6) 



where the baryon number density is ?t-b = "-quark /N c at 
the boundary z = 0. This determines the normalization 
of the baryon number, in the solution ^ of the D7-brane 
system. 

Horizon formation on the flavor D7-brane 

A time-dependent configuration on the D7-brane mod- 
ifies the effective metric which the fluctuations on the D7- 
brane feels. A large field configuration creates a horizon 
of the metric on the D7-brane, which signals the ther- 
malization in AdS/CFT. From the induced metric with 
the background gauge configuration ([5| , here we compute 
the location of an apparent horizon on the D7-brane [28] . 

Let us compute an effective metric which a scalar fluc- 
tuation r\ feels, which is massless scalar meson like pion 
field. By expanding the D7-brane action ([2| for small 
fluctuation 8r\ in the background solution F tz ^ , to the 
quadratic order, we can obtain effective metric g as: 



5' 



dtdzd 3 x i d 3 9 I 



~-^g MN d M 5vdNSv 
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where i — 1,2, 3 is the spatial directions of our 3+1- 
dimensional space-time, and 9 1 (I = 1,2,3) is the an- 
gular variable on the S 3 . (M,N) shows the whole 7+1 
directions on D7-branes, (t,z,i,I). The effective metric 
g can be obtained easily, 
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where G/j is the metric on the unit 3-sphere. 

Given this effective metric, we will now determine the 
apparent horizon, which is defined locally as a surface 
whose area variation vanishes along the null rays which 
is normal to the surface. The surface area at an arbitrary 
point in given (t, z) is 

Surface = J d 3 x t d 3 9 I ^/{U t=1 ^ 3 ~ 9ll ) (II /=1 , 2 , 3 g u ) (9) 
= V/gVol^K^^a - z i R-\2na' 2 )Fl)- 1 ' 2 . 

The (t, z) spacetime has a trivial null vector normal 
to the surface (w', v z ) = (1, —1) since g zz — —gtt, so the 
constancy of the surface area variation along this null ray 
is 



dV s 



surface \dt— — dz 



= o, 



(10) 



which yields 

(8 Z - 8 t ) [z 2 (l - z A R-\2^ol 2 )F 2 z )\ = 0. (11) 

Substituting the gauge field solution (|5| to this, we obtain 
the following equation 



(27ro/) 2 i? 12 - 2z 6 g 2 + 2z 7 gg' = . 



(12) 



If this equation (12 1 admits a solution, it specifies where 



the apparent horizon on the D7-brane is formed. 

Before we proceed, a few comments are in order. First, 
we are calculating the apparent horizon, not the event 
horizon on the flavor brane. Since apparent horizon is 
always inside the event horizon, as long as apparent hori- 
zon never disappears at finite time, we can regard the 
formation of the apparent horizon as a signal of the thcr- 
malization of system. However since the positions of 
the apparent horizon are time-dependent, it is difficult 
to extrapolate the thermodynamical information such as 
temperature, which is determined by the event horizon. 
On the other hand, the apparent horizons are defined lo- 
cally without knowing the late time asymptotics, there- 
fore it has a calculation simplicity. Secondary, we are 
using the Born-Infeld action on the flavor D7-brane to 
determine the effective metric for the various mesonic 
modes. This Born-Infeld form is crucial, since if we have 
used the Yang-Mills form for the D7's degrees of freedom, 
the horizon would have not formed. The reason why we 
need the Born-Infeld form is due to the warping in the 
0, which makes the effective string tension a' finite, so 
it is not appropriate to replace the Born-Infeld action by 
the simple Yang-Mills form. 

Thermalization time-scale order estimation 



From (12), we can order-estimate the thermalization 



time-scale by a dimensional analysis without specifying 
the explicit form of the source function g{t — z). Wc 
would like to consider the chemical potential change 
which mimics the heavy-ion collisions. The function 
g(t — z) at the AdS§ boundary z = is directly related 
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to the time-dependent baryon number according to ([6| 
as g(t) = (2/7r)(27ra') 4 An B (f). 

Suppose that g(t — z) changes like trigonometric func- 
tions from zero to some maximal value during the time- 
scale 1/w. This is like the situation where the chemi- 
cal potential change locally by two baryonically-charged 
heavy ions approaching each other. Setting the maximal 
value of g(t — z) as g ma ,x, we can order-estimate it as 



9(0 



; , g'(£) ~ wg n 



(13) 



Then, a dimensional analysis of (12 1 estimates the ther- 
malization time-scale tth as 



Hh 



(2ira 



M2^>12 \ 1/6 



5max 



1/6 



(14) 



if t t hw < 1 is satisfied. Here being the max- 
imal baryon number density determined as g ma x — 
4(27ra') 4 AnB- If t^w 3> 1 instead, then we obtain, 



tth 



(2W) 2 i? 12 ^ 1/7 



A 



1/7 



(15) 



On the other hand, if has an explicit £ dependence 
like a power-law behavior to approach its maximum, such 
as <?(£) cx £™, with positive n, then 



9(0 ~ 5ma X «)™ , g'(0 ~ ^axW)"" 1 

Again, a dimensional analysis yields, 



^th 



(2W) 2 i? 12 



l/(6+2n) 



A 



l/(6+2n) 



(16) 



(17) 



In summary, in terms of following parameters, inverse 
of the variation timescale of the baryon chemical poten- 
tial w, 't Hooft coupling A, and the maximum baryon 
density «b, the thermalization time scale is written as 



tth 



A 



l/(6+fe) 



(18) 



for given k(> 0), which is determined by how we change 
the baryon number chemical potential. This is one of our 
main results. 

In the following, we present two explicit examples of 
the source function g(£) and show that both examples 
exhibit the generic behavior (18). The first example is 



for a baryonic matter formation, and the second is for 
baryons colliding and passing through each other. 

Example I : Baryonic matter formation 

To understand the time-scale in more detail, let us in- 
vestigate a few explicit examples. The first example pre- 
sented here mimics colliding baryons forming a baryonic 
matter. We start with zero baryon density ne = and 
then increase it linearly in time, for < t < 1/w. It 




FIG. 2: Locations of the apparent horizon in z-t plane writ- 
ten in the unit of 1/w for the parameter is chosen such that 
(A/(27rn|)) 1//6 = 1/w. The blue, red, green curves/lines rep- 
resent the curves/lines (201, t — z + 1/w, and (211, respec- 
tively. 



reaches the maximum at t — 1/w and then it is kept 
constant. This is like two baryonically charged heavy- 
ion approaching each other, followed by a formation of a 
QGP gas with large baryon number. Setting £ = t — z, 
we arrange the function g accordingly as 



9(0 





^max 



W £ 



(C<o) 

(0 < f < 1/w) 
(1/w < 



(19) 



With this choice of the time-dependent baryon number 
density, we can compute the location of the apparent 



horizon from (12 1. The results are 



t = — + - 




and 



2tt 2 / 3 



1/6 



(for t > z + 1/w) . 



The curve ( 20 ) crosses with line t = z + 1/w at z 



where zq satisfies 



1 — z w - 



647r 4 n|zg 



= 0, 



(20) 



(21) 



(22) 



and in the case zqw <C 1, z — zq coincides with the curve 
(21 1. These results are shown in Fig. [2j 

The slight discontinuity between curves (|20| and (|21~ ]) is 



simply due to a cusp of the input source function ( 19 1 at 
£ = 1/w. If we smoothen the source function (19), then 
the two curves (20) and (|2ll) 



are connected smoothly. 
Just for a comparison, in Fig. [3] we also show the loca- 
tion curve for the apparent horizon for a smooth source 
current given by 

1 



9(0 



Smax 2 (1 




tanh(2w;£ - 1)) 

« « 1/2) 
« ~ 1/2) 

« » 1/2) 



(23) 



FIG. 3: A numerical solution for the location of t he a pparent 
horizon for the source current given by g(£) of (23 1, where 
the parameters and the unit are the same as Fig. |2| This is 
a smoothened version of the curve analytically determined in 
Fig-i 



This Fig. [3] is well compared with Fig. [2] 

The emergence of the horizon on the flavor D7-brane 
is seen by the boundary observer through a light prop- 
agation on the D7-brane. Suppose that the point A at 
(tA, 2a) m Fig. [4] gives the earliest delivery of the infor- 
mation of the apparent horizon. Due the fact that the 
light geodesic toward the AdS§ boundary is again along 
the null vector (V, v z ) — (1, — 1), it is clear that the point 
A is determined by the curve ( 20 ) and its tangential out- 
going null line, and it gives 



tA = c A t (X/n 2 B w 2 ) 1/8 , z A = c Az (X/n 2 B w 2 ) L, ° , (24) 

where c At an d c Az are order-one coefficients and given by 
c At = (11V73 - 13) 1 / 2 /(48(2tt) 4 (5V73 - 31) 3 ) 1 / 8 « 0.88 

and c Az = ((5\/73 - 31)/48(2tt) 4 ) 1/8 « 0.33. As a result, 
the thermalization time seen by the boundary observer 
is 

, 1/8 

tth = tA + Z A = (cAt+CAz) 



2\l/8 



A 



1/8 



(25) 



On the other hand, it is also possible that the point 
C in Fig. [5] gives the earliest occasion, depending on the 
values of the parameters in the source function g(£). This 
happens especially if 



A 



n„u> z 



In this case, similarly we can compute the thermalization 
time t t h as 



tth — tc + Zq , 



with 



tc 



1 

w 



zc 



2tt 2 /3 



1/6 



(27) 



(28) 



FIG. 4: Earliest thermalization for the boundary observer oc- 
curs at point A, which is the tangential point between curve 
( 20 1 (blue curve) and the light ray propagating toward bound- 
ary (black line). The boundary observer see the thermaliza- 
tion at point B. We choose the same paramete as Fig. [2] 




FIG. 5: Curves for the parameter chosen (A/(27rn|)) 1//6 = 
(I0) 2 / 3 /w, written in the unit of 1/w. Thermalization occurs 



at point C, which is the crossing point between line ( 21 1 (green 
line) and the line t = z + 1/w (red line). From point C, the 
light ray is propagating toward boundary (black line), and 
the boundary observer see the thermalization at point D. 



where zc is given by the line (21). Therefore 
1 1 



tth = 2z c 



A\ 1/6 1 



B 



Due to the inequality (26 1, 1/w <C 2zc, this yields 

1/6 



t 



th 



B 



(29) 



(30) 



These results (25 1 and (30) are consistent with the or- 



der estimation ( 18 1 in the previous section. 



(26) Example II : Baryons passing through each other 



Let us investigate another explicit example. This sec- 
ond example mimics baryons which first collide each 
other while then pass each other and leave. We start 
with zero baryon density ijb = and then increase it 
linearly in time, for < t < 1/w. It reaches the maxi- 
mum at t = 1/w and then decrease to zero again. This 
is like two baryonically charged heavy-ion approaching 
each other, and then they pass by due to the asymptotic 
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FIG. 6: The curves (20) (blue) and (32) (green) and t = z 
(yellow), t — z + 1/w (red), t = z + 2/w (magenta) and 
light ray (black) toward boundary which is tangential to ( 20 1 
in the unit of 1/w, where the parameter is chosen such that 
(A/(27rn|)) 1/ ^ 6 = 1/w. It is clear for this parameter that the 
earliest apparent horizon is seen at point F, which is similar 
to point B in Fig. [4] 



0.2 0.4 0.6 0.8 



freedom. Setting £ 
accordingly as 



9(0 



we arrange the function g 




(£<o) 

(o < i < i/w) 

■jO (1/w < £ < 2/w) 
(2/w < 



(31) 



where <? max = (2/7r)(27ra') 4 AnB with n-Q is the maxi- 
mum baryon number density. Compared with an ex- 
plicit example I, the location of the apparent horizon for 
< £ < 1/w is again given by the curve (20 1. On the 



other hand, for 1/w < £ < 2/w, the curve is given by 



t = 




A 



z 6 (2Tr) 4 nlw 2 
(for z + 1/w < t < z + 2/w) . 



(32) 



Finally for z + 2/w < t, (12 1 admits no solution, which 
means that there is no apparent horizon in this region. 
In Fig. |6j we plot these curves in the z-t plane. 

Similarly for a comparison, in Fig. [7] we show also 
the location curve for the apparent horizon for a smooth 
source which takes a Gaussian form as 



7(0 = 9, 



exp 



(2wf (i - 1/w) 2 



(33) 



This Fig. [7] is well compared with Fig. [6] 

Similar to the explicit example I, from the curves ( 20 ) 
and ( |32[ ), we can derive the thermalization time tth- From 
(201, wc obtain again (25|). From the other curve (32 1, 



the thermalization point is at the crossing of the green 
and the red lines in Fig. [6j The z value of this crossing 
point, written as zq, is a solution of the equation 



w 



A 



w 2 4(27r) 4 n|ui 



0. 



(34) 



With this, the thermalization time from the curve (32 1 is 
computed as tth = 2zq + 1/w. Therefore, from the two 



FIG. 7: A numerical solution for the location of the apparent 
horizon written in the unit of 1/w for the source current which 
is given by <?(£) of (33), where the parameters are chosen as 
(A/(27rn|)) 1 ^ 6 = 1/w and written in the unit of 1/w. This is 
a smoothened version of the curve analytically determined in 
Fig.© 



curves, we obtain the thermalization time as the earliest 
occasion among these two, 



tth = min < (c At + c Az ) 



A 



1/8 



2z G + - \ . (35) 
w 



Note that in the case where the first term and the third 



term are well balanced in eq. (34), it gives 



ZG 



1/7 



(36) 



On the other hand, if the second term and the third term 



are balanced in eq. ( 34 ) 



1/6 



Therefore, if 



1/w < t th 



(37) 



(38) 



is satisfied, thermalization time-scale ( 35 ) becomes the 
lowest scale 



tth ~ min{ fe=0 ,i,2} 



A 



rikw k 



1/(6+*)' 



(39) 



Again this is consistent with the generic order estimation 
p). 



Comparison with experiments 

We will later discuss the validity and the universality 
of our results in various other theories which are QCD- 
like theories. But before that, it is quite entertaining to 
substitute realistic values of the parameters and compare 
our results with the data, even though our setting is not 
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realistic QCD at this point. For RHIC and LHC heavy 
ion collisions, the baryons are passing through each other, 
so we may approximate them by the example II above. 

First, let us consider RHIC parameters. It is nat- 
ural to assume that tib is twice the standard nuclear 
density tin times the Lorentz contraction factor 7, i.e., 
nB ~ 277JN, where n-^ ~ 0.17[fm]~ 3 . In RHIC experi- 
ments, we have heavy ion Au-Au collisions with A = 197. 
The Lorentz factor is given by the ratio between its en- 
ergy scale E = ^/s NN /2 ~ 100 [GeV] and the mass of 
Au toau, therefore 7 = E/niAu ~ 100. On the other 
hand, at RHIC, the time scale 1/w should be given by 
the time scale of two nuclei passing by through their bod- 
ies, where two nuclei are propagating almost with the 
velocity of light. Therefore, 1/w is well approximated as 
1/w ~ 2^4 1 / 3 /7 [fm/c], where A is the nucleon number 
and 2A 1 / 3, [fm] is the typical nuclear diameter. This gives 
1/w ~ 0.1 [fm/c]. 

With these inputs at hand, we obtain 



Lorentz factor 7 compared with RHIC case in (40 1, (41 1, 



A 
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1/8 
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A 
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1/7 



1/6 



A 2 /*\ 
7 4 < 



A^ 3 X 
7 3 «n 



7 2 «n 



1/8 



1/7 



1/6 



0.24 x A 1/s [fm/c] . (40) 



0.30 x A 1 / 7 [fm/c], (41) 
0.39 x A 1/6 [fm/c] . (42) 



These are all bigger than 1/w ~ 0.1, therefore ap- 
proximately (38) is satisfied. In gauge/gravity duality, 



the 'tHooft coupling A = <?qcd^c is taken to be very 
large. However, since the power of A in tth is small (less 
than one ), tth can not take a large value, even for large 
A. Actually we use A ~ O(10) which is often used in 
gauge/gravity duality for the spectrum comparison. In 
this case, the smallest of these are given by (40), though 



all the scale ( 40 1 , ( 41 ) , and ( 42 ) give the same order time- 
scale. Therefore we obtain the thermalization time-scale 
as 



t th < 1 [Wc] • 



(43) 



It is interesting that this time-scale can be well com- 
pared with the known hydrodynamic simulation require- 
ment t t h < 2 [fm/c] [T3HI7] . We found a rather rapid 
thermalization of mesons. 

Our calculation can also gives a "prediction" for heavy 
ion collisions at LHC. In LHC where Pb-Pb ion collision 
experiments are on-going, the energy scale is bigger than 
RHIC as yswv = 2.7 TeV [2T]- The Lorentz factor is 
7 = E/mpb ~ 1000, due to the center of mass differ- 
ence, and therefore, 7 is 10 times bigger in LHC than 
RHIC. Since A s» 200 is almost the same, these give 
1/w ~ 0.01 [fm/c]. With these at hand, we can compute 
the thermalization time-scale. Due to the difference of 



( 42 ) , the thermalization time-scale is suppressed further- 
more in LHC, and we obtain 



t th < O (0.1) [fm/c] 



(44) 



Again, with A ~ O(10) is used. (Due to the asymptotic 
freedom, Alhc < ArhiCj however the difference between 
Arhic arL d Alhc is tiny therefore we can neglect this ef- 
fect.) This (44) gives a significantly faster thermalization 
time-scale compared RHIC. 

Thermalization for various modes 

Given the calculation of the scalar meson thermaliza- 
tion in the massless Af = 2 supersymmetric QCD, it is 
straightforward to extend the calculation of the thermal- 
ization for other degrees of freedom. 

First, instead of the thermalization of the scalar me- 
son 77, let us consider that of vector mesons Am- A sim- 
ilar computation leads to an effective metric of Am on 
the D7-brane, which gives the equation for the apparent 
horizon as 



(d z - d t ) [z 4 (i - z 4 r 4 (2W 2 )^)] = . 



(45) 



This differs from (11) by just a power in the z factor. 



According to this modification, the thermalization time 
is just 2 1 / 3 times that of the scalar meson. Therefore 
the thermalization time-scale ( 43 ) is almost common in 



order, for various vector meson excitations on the flavor 
branes. 

Next, we consider effects of the quark mass. The quark 
mass m q corresponds to the boundary location of the D7- 
brane, r\(z = 0) = 27ra'm q . This shifts the D7-brane a 
bit. One can compute the full effect of this shift in our 
formalism, but we can give a naive estimate as follows. 
Noticing the fact that 77 comes in the effective metric 
always as a combination (R 2 /z 2 +rj 2 /R 2 ) instead of just 
R 2 /z 2 , our computation presented here for 77 = is valid 
when 



R 2 /z 2 > rj 2 /R 2 . 



(46) 



Substituting the expression of z by zc given by ( [28] ) , we 
obtain 



m q < (^Atib/tt) 1 / 3 . 



(47) 



Even without relying on the large 'tHooft coupling limit, 
this is generally satisfied for light three flavors, for the 
standard nuclear density. Therefore again we expect that 
the thermalization time-scale ( 43 ) is almost common in 



order, even for the various meson excitations with differ- 
ent flavors, such as up, down, and strange flavors. 

Discussion on universality and real-world QCD 

We saw that either in massless or massive Af = 2 super- 
symmetric QCD, the calculations of the thermalization 



time-scales of the various meson modes are always given 
by (18 1. Given this, it is natural to ask to what extent 



our thermalization time-scale ( 18 ) holds for a larger va- 



riety of gauge theories. Since this question is related to a 
possible universality and also to real-world QCD-related 
problems, we shall discuss this question now. 

First, note that even though our setting admits super- 
symmetry, it is also clear that we have never used the 
fermion properties for our thermalization calculations. 
Therefore we expect that our results are not much de- 
pendent on the supersymmetry. 

Next, we shall see that, even with different background 
metrics, our thermalization timescale ( 18 1 is universal. 
We have used the AdS^ metric 0, which represents the 
deconfined phase for the gluon sectors in the conformal 
J\f = 4 theory. However, the conformality of the gluon 
sector in the metric is not important at all since 
our computations of thermalization reply only on the 
asymptotic part [z < zq) of the induced metric, where 
zq is the point where apparent horizon emerge (such as 
za, zc, ze, zq in our explicit examples I and II). Therefore 
we claim that our results hold for other theories where 
IR dynamics of gluons are significantly different from our 
theory. Even if we replace the metric |l]) by some other 
non-conformal metric which does not admit a bulk hori- 
zon, such as a cut-off AdS§ at z > zq, due to the fact 
that our calculations are insensitive to the IR regime of 
the geometry at z > Zq, our conclusion is still valid. In 
this sense, we expect that our thermalization time-scale 
( 18 1 is not only for the J\f = 2 theory, but rather it works 



for a broader category of non-conformal theories. 

The scale zq < O(l) [fm] in the example I and II we 
studied corresponds to the energy scale > 0(200) [MeV]. 
Therefore for any non-conformal theory which admits 
confincment/deconfinement transition for the gluon sec- 
tors at the scale smaller than 200 [MeV], our result (43) 
for the thermalization of mesons is expected to be valid. 
For the real- world QCD, confinement/deconfinement 
transition is ~ 200 [MeV], which is the validity bound 
of our analysis, therefore it is expected that our results 



(43 1 also hold even for realistic QCD. However to confirm 



this, furthermore study is necessary. 

Let us discuss confining gauge theories a bit more. On 
the gravity side, confinement is implemented as a defor- 
mation at the IR region (z ~ oo) of the geometry. If we 
use successful hard-wall models in the bottom-up models 
of holographic QCD, which is with a cut-off of AdSs at 
IR z > zq, then we obtain the same thermalization time- 
scale. One other example is a confining geometry made 
by D3-D(-l) system which has the AdS 5 form at 

the UV region. In the solution the D(-l)'s (D-instantons) 



condense in the bulk and back-react to modify the IR 
r — > region by emitting the dilaton (f> as 



ds 2 = e* /2 x (AdS 5 x S 5 geometry) 



= 1- 



(48) 



Here ds 2 is the metric in the string frame. In the UV 
r — > oo, this reduces to the AdS metric ([I]), as adver- 
tized. The D(-l) charge q (which is proportional to QCD 
instanton charge density) is related to the QCD string 

This breaks the 



-l/r-2 



'qcd- 



tension tqcd, as q/R = ir A 

supersymmetries by half. The universality is valid if this 
factor may not significantly modify the asymptotic 
geometry around z < zq, therefore we need to require 
q < R s z 



With ( 28 1 , this translates to a condition 



TQCD < 



2 5 / 3 A 1 / 6 4 /3 . 



(49) 



Realistic parameters used in this paper show that tqcd 
is at a comparable order with the right hand side, so 
this effect may modify the thermalization time-scale only 
slightly by an O(l) factor. Further study would be inter- 
esting for these. 

Finally we comment on possible generalization of our 
approach for future works. We have conducted the calcu- 
lations with the abelian Born-Infeld action, which treat 
a single flavor brane. In order to treat multi-flavors, we 
need to extend our analysis to non-abelian Born-Infeld 
action. It is interesting to generalize our study to non- 
abelian flavor branes and study the flavor dependence of 
the thermalization time-scale. Since our thermalization 
calculations are insensitive to the IR but are sensitive 
to the UV regime, if we consider totally different bulk 
geometries which do not approach the AdSs metric ([I]), 
our result ( 18 ) is not valid any more. It is quite inter- 



esting to generalize our approach to other theories where 
their UV geometries arc different from ours, such as the 
holographic QCD model by D4-D8 on Witten's geometry 
[2~4l |2"5] , or Lifshitz type of geometries [25] for an applica- 
tion to condensed matter systems. We leave these studies 
for future works. 
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